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ABSTRACT 

We report on the initial analysis of a Herschel / PACS full range spectrum of Neptune, covering the 5 1 -220 yum range with a mean resolving power 
of ~3000, and complemented by a dedicated observation of CH4 at 120 yum. Numerous spectral features due to HD (R(0) and R(l)), H2O, CH4, 
and CO are present, but so far no new species have been found. Our results indicate that (i) Neptune's mean thermal profile is warmer by ~3 K than 
inferred from the Voyager radio-occultation; (ii) the D/H mixing ratio is (4.5±l)xl0"^, confirming the enrichment of Neptune in deuterium over 
the protosolar value (~2.1xl0"^); (iii) the CH4 mixing ratio in the mid stratosphere is (1.5+0.2)xlO"'', and CH4 appears to decrease in the lower 
stratosphere at a rate consistent with local saturation, in agreement with the scenario of CH4 stratospheric injection from Neptune's warm south 
polar region; (iv) the H2O stratospheric column is (2.1±0.5)xl0''' cm"' but its vertical distribution is still to be determined, so the H2O external 
flux remains uncertain by over an order of magnitude; and (v) the CO stratospheric abundance is about twice the tropospheric value, confirming 
the dual origin of CO suspected from ground-based millimeter/submillimeter observations. 

Key words. Planets and satellites: individual: Neptune; Techniques: spectroscopic; Infrared: solar system; Radio lines: solar system) 



1. Introduction 

Neptune's thermal emission has been initally explored from the 
ground in the 8-13 window and in the millimeter range and 
by the Voyager spacecraft in 1989, but detailed views of its 
spectrum had to await sensitive instrumentation onboard ISO 
(see review in Bezard et al. 1999a), Spitzer (Meadows et al. 
2008) and recently AKARI (Fletcher et al. 2010). Altogether, 
these observations have revealed a surprisingly rich composi- 
tion of Neptune's stratosphere, including numerous hydrocar- 
; bons (CH4, C2H2, C2H6, CH3, C2H4, CH3C2H, C4H2), oxygen- 
, bearing species (CO, CO2, and H2O), HCN, as well as deu- 
terium species CH3D and HD. Favorable factors for observing 
minor species in Neptune's atmosphere are (i) its relatively warm 
stratosphere (~140 K at 1 mbar) that enhances IR emission; and 
' (ii) Neptune's large internal heat source that results in rapid con- 
vection updrafting minor disequilibrium species, notably CO, 
up to observable levels. Neptune's submillimeter spectrum long- 
wards of 50 yum has been observed by ISO/LWS (Burgdorf et 
al. 2003), but the signal-to-noise ratio in the data was not high 
enough to reveal spectral features. In this paper, we report the 
first results from observations of Neptune at 51-220 pm (195- 
45 cm"') with the PACS instrument onboard Herschel (Pilbratt 
et al. 2010), performed in the framework of the KP-GT "Water 
and Related Chemistry in the Solar System", also known as 
"Herschel Solar System observations" (Hartogh et al. 2009). 



Table 1. Summary of observations 



Obs. ID 



Start Time 
[UTC] 



T 

[min.] 



Range_^ 

[pm] 



1342186536 


30-Oct-2009 


00:58:36 


116 


51-72*, 


102-145" 


1342186537 


30-Oct-2009 


03:01:48 


133 


51-62', 


150-186" 


1342186538 


30-Oct-2009 


05:22:32 


203 


60-73', 


180-220" 


1342186539 


30-Oct-2009 


08:53:20 


151 


68-85'", 


120-171" 


1342186540 


30-Oct-2009 


11:31:41 


236 


82-102"' 


, 165-220" 


1342186571 


31-Oct-2009 


14:35:00 


82 


118.4 


-120.9" 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



" grating order and filter: k = 2A, 1 = 3 A, m = 2B, n = 1 red 



2. Herschel I PACS observations 

All observations (Table 1) were carried out in chopped-nodded 
PACS range spectroscopy modes ( Poglitsch et al. 2010) l at high 
spectral sampling density. The entire spectral range of PACS 
has been measured at full instrumental resolution AI6A rang- 
ing from 950 to 5500 depending on wavelength and grating or- 
der jPACS Observers Manual 201 Oi l. A summary of the observa- 
tions is given in Table [T] Since blue and red spectrometer data 
are acquired in parallel, several spectral ranges have been ob- 
served in overlap. Given the instrumental spatial pixel size of 
9.4"x9.4", Neptune (2.297" as seen from Herschel) can be con- 
sidered as a point source, and the analyzed spectra therefore 
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Fig. 1. Composite PACS spectrum of Neptune, expressed in line/local continuum ratios. For spectral ranges covered more than once 
(Table [TJ, the observation with the highest resolution has been selected. The region beyond 190 fim is not shown, owing to severe 
mixing of spectrometer orders. The bottom curves are synthetic spectra at appropriate spectral resolution that show the contributions 
of CH4, H2O, CO, and HD lines. 



originate only in the central spatial pixel of the integral field 
spectrometer 



Starting from Level products, the processing of all ob- 
servations was carried by standard PACS pipeline modules 
(Poglitsch et al. 20"T0) l up to Level 1. Individual spectral pixels 
were then scaled onto their common mean in order to improve 
the removal of signal outliers caused by cosmic ray hits. After 
application of an iterative cr-clipping, adapted to the instrumen- 
tal resolution, the remaining data were rebinned onto an over- 
sampled wavelength grid to ensure conservation of spectral res- 
olution. The absolute flux calibration of the instrument and im- 
provements on the relative spectral response function are still in 
progress. Therefore the resulting spectrum was then divided by 
its continuum, to be robust against forthcoming calibration up- 
dates. The composite spectrum is shown in FiglT] It shows emis- 
sion signatures due to CH4, H2O, CO, as well as the R(0) and 
R(l) lines of HD at 112 and 56 /im, seen respectively in absorp- 
tion and emission. At this stage of the data reduction, features 
below ~0.5-l % contrast must be treated with caution. No new 
species are detected at this level. 



A dedicated line-scan high S/N observation of the CH4 
1 19.6 i-im rotational line was also acquired in order to get a high 
precision measure of the CH4 stratospheric abundance. 



3. Analysis and discussion 

3. 1 . Thermal profile and D/H abundance 

Observations were analyzed by means of a standard radiative 
transfer code, in which the outgoing radiance from Neptune was 
integrated over all emission angles. The effective spectral reso- 
lution as a function of wavelength was determined by fitting the 
widths of the H2O lines, whose profile is purely instrumental. We 
initially considered thermal profiles inferred in previous work 
(Marten et al. 2005, Bezard et al. 1998, Fletcher et al. 2010, re- 
spectively from ground-based, ISO, and AKARI observations). 
Below about 0.5 bar, all of them follow the Voyager radio- 
occultation profile (Lindal 1992, see also Moses et al. 2005). 
Above this level, these profiles diverge significantly, showing ex- 
cursions of ~5 K over 10-200 mbar, and even larger dispersion 
(~10-20 K) at lower pressures. Over 50-200 fim, Neptune's con- 
tinuum is formed near the 500 mbar level (T^ ~ 59 K). The HD 
lines typically probe the 10-500 mbar range (peak contribution 
near 2 mbar at line center). Because they show a contrasted ab- 
sorption/emission appearance and because HD is vertically well 
mixed, they provide a sensitive thermometer in this region. For 
HD, we used the same linestrengths as in Feuchtgruber et al. 
(1999). We found that the Fletcher et al. (2010) nominal pro- 
file (their Fig. 5) allowed a much better fit of the HD lines than 
the other two profiles, and achieved optimum fit for tempera- 
tures equal to 0.9xFletcher + O.lxMarten (Fig. 2). This gives 
54.5 K at the tropopause, ~3 K higher than in Lindal (1992). 
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Fig. 2. Neptune's temperature and abundance profiles. CH4 pro- 
files condensing (thick red line) or not (thin red line) in the 
stratosphere are considered. For H2O, profiles A and B are those 
of Feuchtgruberet al. (1997), multiplied by 0.95 and 0.9, respec- 
tively, and the "uniform" profile has a mixing ratio of 0.85 ppb 
above the condensation level. For CO, the profiles of Lellouch 
et al. (2005) and Fletcher et al. (2010) are shown. The black line 
shows the inferred temperature profile. 
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Fig. 3. Neptune's spectrum in the 56.0-56.5 and 111-113 //m 
ranges, showing HD lines at 56.25 yum (R(l)) and 112.1 jum 
(R(0)). Thick red line: model for HD/H2 =9x10"^ and the nom- 
inal thermal profile of Fig. |2] Thin pink lines: same for HD/H2 
= 6x10-^ and 12x10 ^ Green: model forHD/Hi = 9x10"^ and 
Marten's et al. (2005) thermal profile). Note also the water line 
at 56.35 fim, well fitted by profile A in Fig. 2. The upper and 
lower blue lines show models for this H2O profile multiplied 
and divided by 1.5. 



Given Neptune's temperature field as inferred from Voyager 
measurements (Conrath et al. 1998), this is probably related 
to the high latitude (42°S) of the Voyager occultations. Based 
on mid-infrared measurements of ethane, Hammel et al. (2006) 
also found enhanced temperatures (but at sub-mbar levels) com- 
pared to Voyager, a likely consequence of seasonal variability. 
Although the HD lines do not constrain temperatures above the 
1 mbar level (needed in particular for analyzing the H2O lines), 
we retained the 0.9xFletcher + O.lxMarten combination for all 
levels. We determined HD/H2 = (9+2)xl0 ^ i.e. a D/H ratio 
of (4.5+l)xl0"^ (Fig. 3). This is nominally less than but consis- 
tent with the (6.5^j g)xlO"^ value inferred by Feuchtgruberet al. 
(1999) from observations of the R(2) line of HD by ISO/SWS, 
and confirms that Neptune is enriched in deuterium compared 
to the protosolar value (~2.1xl0"^) represented by Jupiter and 
Saturn (Lellouch et al. 2001). We defer a joint analysis of ISO 
and Herschel data to future work. 

3.2. Methane, water, and carbon monoxide abundances and 
profiles 

Methane has been observed in Neptune's stratosphere with 
a range of abundances exceeding the saturation value at the 
tropopause cold trap (e.g. Baines and Hammel, 1994). The PACS 
spectrum shows several rotational lines of CH4 in emission over 
80-160 yum. Thanks to the mild temperature dependence of the 
Planck function in this spectral range, these lines are well suited 
to determination of the CH4 stratospheric abundance. We as- 
sumed a CH4 abundance of 2 % in the deep troposphere, then 
following the saturation law. In the stratosphere, the CH4 pro- 
file was characterized by its high-altitude mixing ratio {qcHt) 
and assumed to follow local saturation below the condensation 
point near 40 mbar Utilizing the Boudon et al. (2010) results on 
the absolute CH4 line strengths and in particular using the high 
S/N dedicated CH4 120 yum line scan (Fig. |4]l, we determined 
qcH^ - (1.5+0.2)xl0"^ consistent with Bezard et al. (1999b) 
((0.5-2) X 10"^) but only marginally with Fletcher et al. (2010) 
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Fig. 4. Methane lines at 119.6 iim (Obs.ID 1342186571) and 
159.3 urn (from Obs.ID 1342186537). Red: Model for strato- 
spheric qcHi - 0.0015 above the stratospheric saturation level 
(thick red line in Fig. 2). Green curves: same, but for qcHt - 
0.0020 (upper curve) and 0.0010 (lower curve). Blue: Model in 
which qcHi - 0.0025 down to ~800 mbar (thin red line in Fig. 
2). 

((0.9+0.3)xlO"^). Because of the progressive increase of the 
continuum level longwards of 100 jum, the CH4 features at 137 
jjm and particularly 159 yum are sensitive to the CH4 amount in 
the lower stratosphere. An alternate assumption would be that 
the CH4 is supersaturated there, as could perhaps result from 
strong convective overshoot. This situation leads, however, to 
unobserved absorption wings at 159 yiim and to inconsistent mix- 
ing ratios for the different lines (Fig.|4]i. A 1.5x10"^' mixing ra- 
tio is ~10 times greater than allowed by the 56 K cold trap, and 
consistent with saturation at 60 K. The most probable origin of 
this elevated stratospheric abundance is that CH4 leaks from the 
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hot (62-66 K at the tropopause) Southern region (Orton et al. 
2007) and is redistributed planetwide by global circulation. A 
combined analysis of the PACS, ISO, Spitzer, and AKARI data 
in terms of stratospheric methane and temperature profile will be 
performed in the future. 

The presence of H2O in giant planet stratospheres, includ- 
ing Neptune's, was established from ISO/SWS 30^5 yum spec- 
tra (Feuchtgruber et al. 1997), demonstrating the existence of an 
external oxygen supply. In Neptune's case, ISO observations de- 
termined a (2^)xl0''* cm"^ column density, but did not estab- 
lish the water vertical profile, a parameter needed to derive the 
rate at which water is removed by vertical mixing and condensa- 
tion and to infer the input flux of water. More than 20 H2O lines, 
encompassing over a range in opacity of more than an order of 
magnitude (~0.2 to 2.5), are detected in the PACS spectrum. If 
uniformly mixed above the condensation level near 1.2 mbar, 
the water mixing ratio is q//,o = (0.85±0.2) ppb, and its col- 
umn density is (2.1±0.5)xl0''* cm"^. Following Feuchtgruber 
et al. (1997), we also considered H2O vertical profiles resulting 
from transport models, characterized by the eddy diffusion co- 
efficient profile (profiles "A" and "B", see Fig. |2|. For a given 
vertical profile, the water amounts we determined from the data 
were identical, to within 10 %, to the values inferred from ISO. 
However, the associated external fluxes vary strongly (1.4x10^ 
cm^^s"' for model A and 9x10^ cm^^s"' for model B). We leave 
the detailed retrieval of Neptune's water profile (including PACS 
targeted observations of several weak lines and a deep 557 GHz 
HIFI observation) for the future. For the time being, an elemen- 
tary analysis based on the integrated linewidths favors profile A 
over the other two water profiles (Fig|5]l, suggesting that the wa- 
ter mixing ratio increases with altitude over 0.1-1 mbar. 
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Fig. 5. Modeled vs. observed H2O line integrated areas for 
the three water profiles of Fig. 2. Line areas are expressed in 
cm ' X % of the local continuum. For each profile, the mean 
rms dispersion (in the same unit) between observed and mod- 
eled areas is given. Profile A provides a better fit to the data than 
do the other two profiles. 

Recent CO observations at millimeter/submillimeter wave- 
lengths (Lellouch et al. 2005, Hesman et al. 2007) point to a 
higher abundance of CO in Neptune's stratosphere than in the 
troposphere. Both studies thus indicate a dual external/internal 
source, with the external source possibly provided by an ancient 
cometary impact. They also provide consistent values of the CO 
tropospheric mixing ratio (0.5-0.6 ppm). However, they differ by 



more than a factor of 2 (1x10^^ and 2.2x10"^, respectively) on 
the stratospheric CO abundance (Fig.|2]l. Support for the Hesman 
et al. value was reported from the detection of CO fluorescence 
at 4.7 yum by AKARI (Fletcher et al. 2010), from which a 2.5 
ppm abundance of CO above the 10-mbar pressure level was in- 
ferred. We find here that the CO lines longward of 150 jum (Fig. 
|6ll instead imply a CO stratospheric abundance of ~1 ppm, in 
agreement with Lellouch et al. (2005). The detailed determina- 
tion of the CO profile will be possible from combined analysis 
of PACS, SPIRE, and new broadband ground-based millimeter 
data. 
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Fig. 6. CO lines at 153-187 /im, compared with models using 
the CO distributions of Lellouch et al. (2005) and Fletcher et al. 
(2010), shown in Fig.El CO fines occur at 154, 163, 174 and 186 
/vm. Other features are due to CH4 and H2O. The bottom curves 
are difference (observed - modeled) plots (shifted by 0.97), fa- 
voring the Lellouch et al. profile. 
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